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Abstract 

We analyze multi-point correlation functions of a tracer in an incompressible flow at scales far 
exceeding the scale L at which fluctuations are generated (quasi-equilibrium domain) and compare 
them with the correlation functions at scales smaller than L (turbulence domain) . We demonstrate 
that the scale invariance can be broken in the equilibrium domain and trace this breakdown to 
the statistical integrals of motion (zero modes) as has been done before for turbulence. Employing 
Kraichnan model of short-correlated velocity we identify the new type of zero modes, which break 
scale invariance and determine an anomalously slow decay of correlations at large scales. 
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When the scale L of an external source of fluctuations far exceeds the scale at which 
fluctuations are dissipated, turbulent cascade appears between those scales. One of the 
most interesting fundamental aspects of turbulence is the existence of anomalies: sym- 
metries remain broken even when symmetry-breaking factors tend to zero. For example, 
time-reversibility and scale invariance of the statistics are not restored even when pumping 
scale goes to infinity and dissipation scale to zero. The mechanism of dissipative anomaly 
(responsible for irreversibility [l|) has been identified by Onsager as due to non-smoothness 
of the velocity field in the inviscid limit P|. It is very much similar to the axial anomaly 
in quantum field theory 3^. Breakdown of scale invariance has been identified recently as 



related to the statistical integrals of motion; in particular, such integrals have been found 
as zero modes of the (multi-point) operator of turbulent diffusion in the Kraichnan model 
of short-correlated velocity P,0,laSyi- 

Here we consider passive scalar in a random incompressible flow and ask whether the 
symmetries are restored at the scales far exceeding the pumping scale. It is straightforward 
to see that time reversibility holds (to put it simply, everythingwhich is pumped is dissipated 
at smaller scales and there is no flux towards larger scales) [S*. At the level of the second 
moment or spectral density of the scalar, an equipartition takes place at wave-numbers 
smaller than like in other systems with a direct cascade 0, ^. Here we flnd four-point 
and higher moments at large scales and discover that the scale invariance may be broken. We 
employ the Kraichnan model and demonstrate that the breakdown is also due to statistical 
integrals of motion (even though very much different from the zero modes breaking scale 
invariance in turbulence). We thus flnd a possible link between an anomalous scaling in 
equilibrium systems (e.g. in critical phenomena) and in turbulence. 

Consider the passive scalar 6{r,t) carried by the velocity v(r, t) and pumped by 0(r,t): 

dtO + {vV)e = (f) . (1) 

The characteristics of (0) are called Lagrangian trajectories and deflned by (9tR(r, t) = 
v(R, t) and R(r, 0) = r. Integrating ([T)) by characteristics we get 9{r, 0) = ^^(R(r, — t), —t) + 



0(R(r, ti),ti) dti. Correlation functions, Fn = (6i . . . 6nj, are obtained by averaging over 
velocity and pumping denoted by brackets and over-bar respectively: 



F^={9 (Ri(-t), -t)...e (R„(-t), -t))+ Jdh.. .J rft„(0(Ri(ti),ti)...0(R„(t„),tn)) .(2) 



Here Rj(t) = R(rj,t) and 6i = 0(rj,O). We assume that the pumping is white Gaussian 



with a zero mean and the variance 0(ri, ti)0(r2, ^2) = x(ri2)5(^2 — ^i), = — rj. 

One can express the correlation functions via the muhi-particle propagators. For example, 
assuming zero conditions at the distant past and space homogeneity, one gets 



Foir) 



° dt f P(n,r,t)x(K)dR 

-00 J 



(3) 



where P(R, r,t) is the probability density function (pdf) of Ri2(t) provided Ri2(0) = r, 
Riiit) = Ri(t) - R,(t). 

We use the Kraichnan model [a, where velocity is Gaussian with the zero mean and 
the variance 



(fa(ri, t)f/3(r2, 0)) = 6{t) [Ko6a(^ - Kal3{ri2 



(4) 
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Translation-invariant propagators satisfy closed differential equations in the Kraichnan 
model 0,0: 



(dt + Cn)P{R, r, t) = 0, P(R, r, 0) = 6{R - r) 



(5) 



where R = (Ri . . .R„) and 2£„ = I] -^a/3(i"ji) ViaVj/3. In particular, the pdf of the two- 
particle separation [-P(R, r, t) integrated over angles] is expressed via the modified Bessel 
function (r^o = d/^y — 1) 

Vo{R, r,t)= ^ ' exp 
D-f\t\ 
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giving the Richardson law of separation W ~ D\t\ 10|. 

We consider 7 > (the case 7 = has been described in Q]), which corresponds to 
Holder exponent of the velocity smaller than unity j^. In this case, the characteristics are 
non-unique which is seen, for instance, from P(R, 0, t) 7^ 5(R). Respective loss of information 
leads to scalar dissipation, which balances pumping by d) and provides for a statistical steady 

n 

state of the scalar 8]. Spatially non-smooth velocities are produced by fluid turbulence 
so that all scales considered in this paper are assumed to be less than the outer scale of 
turbulence that is the correlation scale of velocity. The equation (jS)) implies corresponding 
equation on 

£„F„= x(r.,)Pn-2({rfc,}) . (6) 



The simplest steady state described by (jH)) is the thermal equilibrium with Gaussian proba- 
bility density functional oc exp [— / d'^dr/ (2T)]. Such state has F2{r) = T6{r) which requires 
also zero correlation length L of (p: 

TK^p{v)V^Vp5{y) = -x(r) . (7) 

That is if the limit L — > is taken in such a way that F2 T6{r), the scalar statistics 
becomes Gaussian and scale-invariant. The Gaussian anzatz, F^ = J2 F2{rij)F^_2/n solves 

m, 

in two cases: i) F2{r) = T6{r) (in the compressible case as well) and ii) 7 = 2 where Kaj3{r) 
is r— independent. 

Of course, delta-function is an idealization so let us open a Pandora box of anomalies 
by allowing (p to have a finite correlation scale L. We assume that at r ^ L the function 
x(r) decays faster than any power (say, exponentially). We consider velocity field correlated 
at different points (i.e. 7 < 2). Then, the scalar statistics is no longer Gaussian and very 
different at the scales smaller and larger than L. 

Let us first remind the properties of the turbulent state realized at r < L. A nonzero 
flux x(0) makes the pair correlation function non-smooth at zero: F{0) — F{r) oc r'^xi^)- 
The scahng properties are characterized by the structure functions Spjr) = {[9{r) — ^^(0)]") 
which are given by the zero modes of the operator Cn I13I I. Zero modes are 

functions Z{Ilij{t)} of particles coordinates conserved on average ja ll^- The structure 
functions are determined by the so-called irreducible zero modes (that involve coordinates 
of all the points), their exponents being din Sn/d\nr = n'y/2 — A„ with the anomalous 
exponents given by A„ = n(n — 2) (2 — 7)/2(c/ + 2) in the perturbative domain, A„ ^ nj. 



One can express all the zero modes via the distance, R = y + . . . + -Ri„, and the angles, 
Cj = Rii/R, in d{n — l)-dimensional space: Z{Rij} = R'^f{e). The conservation, CZ = 0, 
in a Lagrangian language means that the growth of the radial factor is compensated by the 
decay of the angular function. The zero modes of the hermitian scale-invariant operator 
come in pairs satisfying the duality relation: for every Z„ = i?'^"/„(e) there exists = 
^7-d{n-i)-cr„ j^^^-j There is the difference between the modes with positive as (that 

contribute at small scales) and with negative as (expected at large scales): CnZ'n produces 
delta-function at the origin or its derivatives (contact terms) rather than zero. A physical 



reason is a nonzero probability of initially distant particles to come to the same point 
For example, C2R'^~'^ oc 5(R) at 7 = 2. 

After discovering that the zero modes Z„ enter the solution of (jH)) at small scales, it is 
straightforward to ask what zero modes contribute at large scales where I/„-F„ = 0. Natural 
candidates to contribute F„ at r > L are Z'^ which decay with the distances. And indeed 
this is true for F2 for which the following representation holds 



D{d--f) Jr D{d--f) 

Provided qo = J x{^) dr 7^ 0, zero mode r'^~'^ dual to the constant dominates F2 at r > L. 
Such duality between r > L and r < L does not work for higher correlation functions. 
Considering, for instance, F4 one observes that the source appearing in © is a linear combi- 
nation of functions independent of one of three on which F4 depends. Since the action of 
Cn on Z'j^ produces contact terms vanishing where any is non-zero, then Z'^ do not satisfy 
the correct boundary conditions at r ~ L and thus cannot contribute F„ at large scales. One 
can show that Z'„ occur in the solution only if irreducible terms are present in the pumping 
correlation functions That means that a Gaussian pumping provides for Z'^ only at 

small scales. Here we show that the reducible terms in the pumping correlation functions 
spawn a new type of zero modes at large scales, which correspond to an axial (rather than 
spherical) symmetry in R-space and to the different duality relation, a'^ = — nd/2 — ak- 

Nonzero go means infinite / F2{r) dr i.e an overall heating. We also call go charge to 
invoke an analogy with electrostatics (literal at 7 = 2 when £ is a Laplacian). When go 7^ 0, 
F2{r) decays by a power law at r ^ L and high-order functions obey normal scaling though 
statistics is non-Gaussian We consider now a pumping with go = 0, like in (|7j), when 

F2{r) decays faster than any power. On a particle language, the distance R(t) in Q explores 
the sphere i? < L in a completely isotropic manner with the result being proportional to 
/ x(r)c/r. Such a symmetry cancellation makes pair-correlation small at large distances. We 
shall show now that mutual correlation between different pairs of distances makes for the 
power-law decays of multi-point correlation functions (interpreted as dipole and quadrupole 
contributions) . 

The identity Q allows us to obtain two important conclusions when all Vij ^ L. Consider 
first t <^ [minrij]'^ /D, then the last term can be neglected since the correlation function of 
(J) is exponentially small. Averaging over velocity at times earlier than —t of which Rj(— t) 



are independent by zero correlation time of v, one derives that the correlation functions at 
large scales are statistical invariants of Lagrangian dynamics: 

F„(ri, . . . r„) ^ F„(Ri(-t), . . . R„(-t)) (8) 

This is equivalent to the above statement that F„ is a zero mode of Ln at large distances. 
On the other hand, the steady state values of the correlation functions can be obtained by 
taking in the limit t — oo when the first term vanishes while the second one saturates. 
For Gaussian statistics of this can only occur due to the rare events where pairs of particles 
traced backwards in time come within the distance L from each other. The distance between 
the pairs remains large with overwhelming probability. Now one can choose intermediate 
times such that (jH)) still holds but Rj(— t) are already such that the inter-pair distances is 
much smaller than the distance between two pairs. This suggests that the main features 
of, say, F4 can be captured by considering the special geometry of two distant pairs. In 
particular, this geometry will determine the scaling of a coarse-grained field, see ()16|) below. 
The overall scaling cr„ of can also be inferred from this particular case 

We consider F4(x, y,z), where x = ri2, y = and z = (ri3 + r24)/2 at x, y <C z. It 
is expressed via the pdf P(X, Y, Z, x, y, z, t) of X(t) = Ri2(t), Y(t) = R34(t) and Z(t) = 
(Ri3(t) + R24(t))/2 conditioned at X, Y, Z(t = 0) = x, y, z: 

/•oo r 

F^^Jdt J P{X,Y,Z,^,y,z,t)H{X,Y)dXdYdZ. (9) 

Here H{X,Y) = xiX)F2iY) + xiY)F2iX) decays where max[X,Y] > L because go = 0. 
The integral is determined by t ~ [max{x,y, L}]'~' /D ^ z^/D, that is the distance Z does 
not vary much. For such times, one can assume X <^ Z, Y <^ Z so that X and Y evolve 
approximately independently. The lowest order approximation of the propagator is via two- 
particles pdfs, 

J dZP{X, Y, Z, X, y, z, t) ^ P(X, X, t)P{Y, y, t) = Po , 

and it produces the reducible part of F4 that is F2{x)F2{y), which is exponentially small 
if max{a;, > L. The next-order corrections decay as powers of x/z and y/z and thus 
dominate the correlation function if x > L or y > L. One can find that correction writing 
in © £4 = ^ + ^1, ^ = K^/3(X)Vxyx^ + K^p(Y)VYyY0, see Appendix El In the first 
order in Li one derives 

Pi = /* dt'dx.'dy'dz'P(X, x', t')P(Y, y', t')S(Z - z')LiP(x', x, t)P(y', y, t)S(z' - z) . 

JO 
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We now plug that into ©, use x = -^2-^2, integrate by parts and obtain the first-order term 

in (x/z)^, {y/zY': 

/•oo _ - 

Fl = dtd^'dy'x'Jf,F^{x')F^{y')P{^',^,t')P{y',y,t') [x%d''K^p{z)/ dz.dzs 
The integration over the directions of x' and y' gives V2 which is P(R, r, t) integrated over 

R 

angles with the second angular harmonics. It was found in [6] that V2 has the same form as 



Vo with Tjo replaced by 772 = 7 ^\J{,d — '^y + Sd{d+ 1 — 'y)/{d — 1). Time integration also 
can be done explicitly and we obtain (see Appendix |X1 for the details of calculations) 

L D ^''^■'y^y' dz dz y^^'^2(^') Jdy'F^iy'){x'yr^'+'^'Q,,-i/2{w), (10) 

where Qri2-i/2 is the Legendre function of the second kind and 8w = 
{xyx'y')^'^/'^[{x^ + y'^ + x'^ + y'^Y ~ ^{xx')'^ — A{yy')'^]. The function Fl approximates 
F4 — F2{x)F2{y) in the whole region z ^ L, x,y <^ z irrespectively of L. At x,y ^ L the 
correlation function given by a zero mode of £4 reduces to the zero mode of Lq: 

pi {xyYql .... d'K^pj^) 

^ (^. + ^.)(4.+2.-.)/.^"^7y/^^^^;;^' (11) 

where 26 = 7772 + 7 — and q2 = J x'^'^^xi^)^^- We thus have shown that, despite the fast 
decay of F2 at large scales, F4 decays as a power-law there and found the scaling exponent: 
a4 = 2d + 26. The dual zero mode of Lq with the scaling exponent a'_^ = •j — 2d — (first 
found in 1^) appears at x,?/ -C L: 

Z oc {xy) XaX^ypys ■ 

One can easily verify for any vector v that x^~'^[d{'x ■ v)^ — x'^v^] is an anisotropic zero 
mode of the Lagrangian evolution operator in the lowest order in the distance x between 
two particles. The above derivation makes explicit the physical origin of the anomalous 
scaling thus found: The main contribution into F4 at large scales comes from Lagrangian 
trajectories that were in the past separated into two distant pairs. If the other pair would not 
be there at all, the average pair prehistory would be completely isotropic that is determined 
by Q'o = / x{^)d^ which vanishes. However, the presence of the other pair leads to a preferred 
direction in space allowing each pair to exploit the quadrupole term q2 = J x(x)x^"'''^(ix. Since 
the correlations between the pairs decay only as a power law, the quadrupole contribution 
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dominates the correlation function. The possibihty for more than two particles to explore 
non-isotropic configurations (which show themselves via non-isotropic zero mode) exists for 
higher correlation functions as well. To give an example of power-law correlations in high 
moments, let us derive the explicit expression for n— point correlation function in the limit 
of small ^ = 2 — 7. Besides reproducing the correct scaling it shows how non-trivial the 
angular structure of the correlation functions is at large scales. Since the statistics of is 
Gaussian at ^ = we analyze the cumulant = lim^^or2n/^ (see Appendix |B|): 



EV^,rL = F(n)Evr,V: 

k {ij} 



k=l 



(12) 



where the sum is over all permutations of indices, F{n) is an n— dependent constant and 



Ja/3(r) = Din r6ai3 - r^rpr "^D/ {d - 1) 



(13) 



The solution at rjj ^ L is as follows 



2ri ~ 



-1) 



2„n-2 



71+1 



AT{nd/2 + A)F{n)qlql 
■''d/W-\d - l){d + 2)^{2dy 



(14) 



{ij} 



r ■ 

*2 ' ■ ■' «n 



c 



,nd+4: 



4.} 



ft 

k=l 



'l'2k-\ 



'«2fc 



where qi = J x'^xi^) and a homogeneous function of zero degree /(ri, .., r2n) is expressed 
via the hypergeometric functions denoted Fijk = F{nd/2 + 4 + i, j + 1, d/2 + 2 + k,u) of 
the argument u = —2z^/c^ (c^ = C|]^ 2 2n}) ^"^^ scalar products ai = (ri2 ■ r34) /ri2r34, 
02 = (ri2 ■ z)/ri2^, as = (r34 ■ z)/ruz: 

u{nd+i 



/vi2 ^2 [ 
r_ '12' 34 p I S^2^2 
/ — ) -^000 + 0020-3 



d + A 



121 • 



Foio + F(^ 



000 



af [2Foio + {d- 2)F( 



000 



+2 (^al + afj (Foio - Fooo) -2010203 



Au{nd+S] 
d + A 



Fi2i + (rf-4)(F( 



010- 



000 J 



(15) 



One can check that in the limit of large z ()15p reproduces ()11|) at ^ (where 6-^2). Let 
us stress that the cumulants r„ are nonzero only when the balance equation ((7j) does not 
hold. If to replace 6{r) in (|7j) by a function with a finite L then go remains zero but dipole 
and quadrupole moments, qi and q2, are generally nonzero an d (TTTl4ll5j) appear. The proof 
that (T„ = nd/2 + 26 for finite ^ will be published elsewhere |l5|. 
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To study how statistics changes with the scale one defines the field coarse-grained over 
the scale r: 9i = J^,^^9(r') dv' /i'^^'^. A fast decay of the pair correlation function together 
with power-law correlations of higher functions, (|llll4ll5p . mean that the statistics of 6^ 
approaches Gaussian anomalously slow as £/L oo. Indeed, while the second moment {6j) 
tends to a constant (exponentially) fast, the fourth cumulant may decay slower than 
(that one would have for fields with a finite correlation radius): 

{{ej)) oc c^''~^dc J Fiz'^-Hz oc {L/if\ (16) 



Here, the reg ion c ~ L ^ 2; ~ £ gives the main contribution if 27 < (see [la] for details 



and also 



ISl]). Similarly one can derive ((^^(0)6'^(r))) oc r "^^ . Note that Fl gives zero 



contribution and one must account for Fl oc z~'^"' . 

Since o"„/n depends on n then the statistics is not scale-invariant at r > L when go = 0. 
If go 7^ then the zero modes found here provide for an anomalous scaling of sub-leading 
corrections, more similar to what is generally observed in critical phenomena. What we 
believe is of importance here is that we trace this anomalous scaling to zero modes. That 
shows that at least in passive scalar problem, the mechanism of an anomalous scaling is 
common for turbulence and thermal equilibrium and raises an intriguing possibility that in 
other problems in quantum field theory and statistical physics one can relate anomalous 
scaling to statistical integrals of motion. 
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their results prior to publication. This research has been supported by the Minerva grant 
8464 and the EU Network " Fluid mechanical stirring and mixing: the lagrangian approach" . 
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APPENDIX A: CALCULATION OF F4 AT A SPECIAL GEOMETRY 

The expression for the four-point correlation function -F'4(ri, r2, rs, r4) follows from Q 

X[Rl2{tlMR34{t2)] + 





dtidt2{ 

-00 



X[Rl3ih)]x[R2Ait2)] + x[RuitiMR23it2) 



(Al) 



Here we consider calculation of F4 for the case where the absolute values of the inter-pair 
distances x = ri — r2 and y = r^ — r^ are much smaller than the absolute value of the distance 
z = [ri + r2 — r3 — r4]/2 between the centers of mass of the pairs. We shall also assume that 
z ^ L. It is easy to verify that the first summand dominates ()A1|) in this case. Indeed, the 
other two terms are contributed by the events on which z shrinks to the distance of order L. 
These events have smaller probability than similar events involving x and y justifying the 
above statement. Finally, using zero-correlation time and time-reversibility of Kraichnan 
model and employing H{x,y) = x{^)F2{y) + x{y)F2{x) introduced in the main text we 
obtain 



/■oo nco 

/ dt{H{R,2it),Rum= dtH{X,Y) 
JO Jo 
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xP(X, Y, Z, X, y, z, t)dXdYdZ, 



(A2) 



where we introduced the joint probabihty density function P(X, Y, Z, x, y, z, t) of X(t) = 
Ri(t)-R2(t), Y(t) =R3(t)-R4(t) andZ = [Ri(t) + R2(t)-R3(t)-R4(t)]/2. In Kraichnan 
model this function satisfies a closed evolution equation 



dtP=[Lo + L^]P, P(t = 0) = (5(X-x)5(Y-y 



dX^dX, 



Li 



1 



Kap I^Z + 
_ X + Y' 



X- Y 



+ K, 



X + Y\ 



a/3 



+ K^p I^Z + 

+ 



Y-X 



dZ^dZ, 



K, 



0/3 



Z + 



X + Y 



-K^p Z 



-K^p ( Z + 



X + Y 
Y-X 
X + Y^ 
Y-X 
X + Y^ 
X- Y 



K. 



aP 



Z + 



X- Y 



dX^dYp 



K, 



X + Y 



aP 



92 



dX^dZp 
+ K^p (z 



X- Y 

2 , 

1 ^ 
2 

Y-X 



Z + 



X + Y 



dYadZp 

Next we observe that the integral is determined by times of the order [maxx, y, L]^/Z? 
(we assume go = here so that F2 is small outside L). During these times which are much 
smaller z'^ / D the distance Z does not vary much while X and Y remain much smaller than 
Z. The main dynamics at these times is just independent evolution of "fast" degrees of 
freedom X, Y governed by the operator Lq. Considering Li as a perturbation we construct 
P as a series 



n=0 



dPn 

dt 



LoPn + LiPn- 



(A3) 
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where Pq is given by Po(X, x, t)Po(Y, y, t)5(Z — z). Here Po(X, x, t) satisfies the initial 
condition Po(^ = 0) = 5 (X — x) and the evolution equation 

(A4) 



aPo(X,x,t) 92Po 



9t '9X„9X^ 
The explicit solution for P„ with n > reads 

P„ = /* dt'dx.'dy'dz'Po{X, x', t - t')^o(Y, y', t - t') 

JO 



x5(Z-z') Li , , ,P„_i(x',y',z',t') 

V / x,y',z' 

Above the subscript of Hermitian operator Li signifies that it acts on variables x', y', z'. The 
asymptotic series for P integrated according to produces asymptotic series P4 = Y^F"', 
where _F"+i/_p" tends to zero as 2 — * 00. To write this series most concisely we note that 



J dXdYPo(X,x',t')^o(Y,y',t' 



x{X)F2{Y) + x{Y) 



XP2(X) 



J rfXrfYPo(X,x',t')^o(Y,y',t' 



F2{Y) 



dX^dXp 



dYadYpi 



-— I rfXt/YPo(X,x',t')i^o(Y,y',t')^2(X)P2(l^), 



where we used the equation satisfied by F2. Using the above we find the following expression 
for P" 



F"" = j dtj P„(X, Y, Z, X, y, z, t)H(X, Y)dXdYdZ 

dtd^'dy'dz'F,{x')F2{y')(L^) , , , P„_i(x', y', z', t) 



X' ,y' ,z' 



Using the expression for P„ we find the expression for P" with ri > 2 in terms of P^: 



P'^= / dtdx.'dy'dz'F^LiPn-2{x.',y',z',t). 
Jo 

In particular, for P^ we obtain 

p2= / dtdxV^o(x,x',0^o(y,y',0^i^'(x',y',z) 

Jo 

We pass to study P^ which can be written as 

F' = J dtd^'dy'Poi^, x', t')Po{y, y', t')x'J'^F^{x') 



(A5) 



(A6) 



xP^'d/' 



Z + + Z -— 



-K„/3 z H — - Ko,f3 z H — 



(A7) 
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Expanding the difference in brackets at large z we obtain 



dtd:>i'dy'Po{yi, x', t')Po{y, y', t')xJ'pF2{x') 



dz.ydzs 



+ z 



-7-2 



It is convenient to perform angular integration by substituting the integrand by its value 
averaged independently over the directions of x' and y'. We introduce functions averaged 
over the directions of x': 



(^X ^X^Pq[x, X , t^) angle S^-yA -\- X^X-yB . 

Taking the trace of the above equation and multiplying it with XaX.y leads to 

(Po(x',x,t))„„gie = dA + B, 

{{x' ■ xf Po(x', X, t))angle = A + B. 



(A8) 



(A9) 



Isotropy of velocity statistics implies that the LHS of the above equations are independent 
of X so that A = A{x,x',t), B = B{x,x',t). It follows from the above relations that 

(1 - d)B = { 1- d{x ■ xf Po(x, x', t)) angle- 



(AlO) 



Using incompressibility we can write the answer in terms of the function B solely 



dz^idzs Jo 
B{y,y',t)x'y'Fi{x')F^{y'). 



dtd'x'dy'B{x, x', t) 



(All) 



The function Bix^x^t) can be found by expanding Po(x, x',t) in Jacobi polynomials in 
X- x'/ (xx'), see j6|,|l5|. One finds that B is the coefficient near the Jacobi polynomial of the 
second-order and it is given by 



P(r, r\t) 



Sd-W{d-l)^t 



exp 
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'y^{d-l)tD 



j^{d-l)tD^ 



d ^2^8^(^+1-7) 



7^(^-1) 



Up to a proportionality factor the above coincides with Pi from 6]. The difference due to 



the misprint in [6| is discussed in more detail in [15|. The integral over time is given by 

{xyx'y')^/^-'^/^ 



B{x', X, t)B{y', y, t)dt 



(27r)(^'^-i)2(rf-l)P> 



Q 



V2-1/2 



{x^ +y^ +x'^ +y'^y -4:{xx'y -4:{yyy^ 



]{xyx'y')^/^ 
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The above leads to the result pOj) from the main text. The analysis of the expression 
proceeds along the following lines. We first consider where x ^ L and y ^ L. In this 
situation one can assume x' x and y' ^ y and use 

\x^ + y^ + x'^ + y'^f - A{xxy - A{yyy\ 



y ^{xyx'y')'^/'^ 
r(r/2 + l/2)r(l/2) (^{xyx'yy'^Y^^'^ 



(A12) 



2m+i/2Y{r]2 + l) \ {x^ + y^f ) 

where we used the large argument asymptotic form of the Legendre function of the second 
kind. Substituting the above into the expression for we find 

^1 _ {xyYxgX^y/sys d^K^/sjz) 2^''^ 

In the opposite limit x <^ L and y <^ L one finds zero mode dual to the above. Indeed, in 
this situation one can use the inverse inequalities x <^ x\ y <^y' and 



■ (x^ + y^ + x'^ + y'^f - A{xxy - A{yyy 



One obtains 



V S{xyx'yyl^ 
r(r^2 + l/2)r(l/2) ( ^{xyx'yyiy ^'" 

2'?2+i/2r(r/2 + 1) V (^'^ + y'^f J 



^ 2^-'-r(,, + i/2)r(i/2) y^y fK^^i-) 



(A14) 



X 



/ dx' \ dy'F^{x')F^{y') 
Jo 



{x'y'Y 



Note that now the dependencies on x' and y' in the integral cannot be separated. 

It is easy to see that F^ vanishes ii x oi y or both are zero. Indeed, if x = then by 
isotropy -Po(x', 0, t) = Po{x', t). The angular integration then produces zero by identities like 

Xo^K^^ (z + dS = 0. (A15) 



lx=R \ Z 

The same argument shows that F^ averaged over the directions of x or y or both vanishes 
as well. As a result such objects as ((^^(O)^^(r))) and {{9^)) are determined by F^. These 
are calculated in ll 
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APPENDIX B: n-TH ORDER CORRELATION FUNCTIONS IN THE LIMIT 
OF SMALL ^ 

We remind that passive scalar statistics is Gaussian at ^ = 2 — 7 = 0. Thus it is 
convenient to discuss the hmit of small ^ in terms of irreducible correlation functions that 
are proportional to C, in this limit. 

The analysis proceeds as follows. We first derive equations satisfied by the irreducible 
correlation functions r„ in the limit of small ^. We show that r„ solves the potential 
problem (electrostatics) in the space of nd dimensions with the source concentrated within 
a generalized cylinder. After the integration over the directions parallel to the cylinder axis 
we express r„ as an integral of a kernel with locahzed source. The large-scale asymptotic 
expansion of r„ is then derived exactly in the same way as the multi-pole expansion in 
electrostatics with the coefficients given by two quadrupole and n/2 — 2 dipole moments of 
the pumping correlation function. We proceed to the calculation. 

1. Equations on r„ in the limit of small ^ 

We stressed above that the statistics of 9 is nearly Gaussian at small ^ so that r„ are 
a small correction to the reducible part of the correlation function which is expressible via 

F2. Thus the source appearing in the equation on the leading order term in r„ should be 
expressible in terms of F2. It is simplest to show this for r4 = F^ — Ff = F4— F2{ri2) F2{rs4) — 
F2{ri3)F2{r24) — -^2(?^i4)-^2(?^23)- Usiug the equation on F4 it is easy to show that 



where the sum over ijkl signifies the sum over all 24 permutations of 1234. The boundary 
conditions on r4 are vanishing with the difference of any of its arguments |rj — rj\ tending 
to infinity. Equations on the irreducible correlation functions of order higher than 4 involve 
lower order correlation functions. Let us consider the equation on Fq defined by 



i>j ijkl - 



Kap{rk - rj) ViaVji3F2{vi - Yk)F2{vj - ri), 



(Bl) 



(B2) 
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Using the equation on r4 one finds for Fg 

1 



i>j,k^i,j 



-Ka/3{rj - Tk 



(B3) 



One observes that the term that involves r4 is of the second order in non-Gaussianity that 
is at small ^ it vanishes as Thus this term is negligible to linear order in ^ and in this 
order pair-correlation function fully determines the source in the equation. This holds for 
r„ with n > 6 as well. Indeed, to linear order in ^ one can neglect all products of F— s that 
enter the definition of r„ (say, in the definition of Fg one can neglect F4 x F4 term. To linear 
order in ^ the definition of F„ reads 



n/2- 



^ „ ^2k ({l"ii_2,..,2fc} 



k=2 {ij} 



{2k)\(n-2k)\ 



X 



^n-2k ({^i2fc+i,..,n}) ' 



(B4) 



where Yj{ij} runs over all permutations of 1,2, ..,2n and the factor [(2A;)!(2n — 2A;)!]~^ lifts 
multiple counting of terms. Using the equation on it is easy to see that the equation on 
to linear order in ^ coincides with 

n-l 



^271^2x1 + 



n 



E 



fc/3 



i>j,k^i,j 



K. 



a/3 1 



n— 1 



F2(r,-r,)F2^_2({rfc^,,,}) 
£fe}(2^)K2n-2;t)!^''^'' ({r.M....^.}) 

^2n-2k ({''j2fc+i,..,2n}) — 0- 

Thus the equation satisfied by F2n can be written as follows 



(B5) 



+ C{n) E 



n F2 - r,,. 



k=l 



0. 



The value of the constant C(n) is not essential for our purposes. To derive the equation in 
the first order in ^ term F^^ = lim^^o ^2n/C we substitute C2n by its value at ^ = 0. Using 
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the fact that I]i>j VjaVja = — I]Vi/2 when acting on translation-invariant functions we 
find 



n 

n H (ri,,_, - r, 



fc=i 



(B6) 



where F{n) is an n-dependent constant, H = F2{$, = 0) and Jais = c^?-K^a/3(r, Ok=o- Since 
constant part of J^fs cancels in flB6|) we will use the equivalent expression 



J,^(r) = D(rf-l)lnr5,^-D^. 



(B7) 



2. Solution of the equation on and its large-scale asymptotic expansion 



The solution can be expressed with the help of an auxiliary tensor Tap defined by 

2n 

V^Ta^(ri, r2, .., r2n, Xi, .., x^) = Jaisi^i — r^) 



One has 



Using 



n=l 

n 

X n (r2fc-i - r2fc - Xfc) . 

k=l 



/n n 
n ^(x,) n 

itji k=l k=l 

Taf3{^ii ) Tjj , Tjg , Fj^j , Tjg . . , r^j^ , Xi , X2, X3, . ., Xfc) 
+^a/3(l"i25 ^'44? ''js' ^i5--y ^hnJ -^l; -^2) X3, .., Xfc) 
Ta/B 1 ^ii 1 ^13 1 '^141 '^is ■ ■ 5 ^i2n 1 ■'^l 5 -^^2) X3, .. , X/j^ 



,2 T{d/2-l] 



4v]-'i/2|r 



^Hd-2 



[B8) 



(B9) 
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we find 



i\2 



nd—1 



X 



{nd- 2)1 



2n 



Ef (r. - r: 
r((n + 2)rf/2- 1) 

22+(n-2)(i/2y^(n+2)cl/2 



i = l ^ no 

Jaf^{r[ - rg)^ (r'l - r'^ - Xi) S (r^ - - X2) 



(n+2)cZ/2-l ' 



where we used 



r(n/2 - 1) 



J [Er=i(x.-x02]"/2-i 

r((n-m)/2- 1) 

(n-m)/2-l ' 



e: 



'\2 



which is readily proved as a relation between Green's functions of a laplacian in different 
dimensions. Introducing auxiliary vectors zi = ri, Z2 = r2+xi, Z3 = rs— r, Z4 = r4+X2— r, 
where r = r'^^ — r'^ we have 



-T, 



_ r((n + 2)d/2-l) 

a/3 - 22+(n-2)<i/27^(n+2)(i/2 J ^^'^^PK'^) I """l 



(2r;-Ez./2)%E-.-(E-0'/4 

- l-(n+2)(i/2 



A;=3 



(BIO) 



The integral over r'^ is easily calculated and we find 

^ r((n + l)ci/2 - 1) 

~ 2^-d/2^{n+l)d/2 



X 



:r - b)' + 



l-(n+l)d/2 



(Bll) 



where we introduced b + r4 + X2 — ri — r2 — Xi and = Efc=i (i'2fc-i — r2A; — x^)^. It 
is convenient to introduce 



Fap{h.,(? lO) = r(a) j dr- 



(r - b)' + c2 



(B12) 
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In terms of this tensor the answer takes the form 



^2n — 5Z 



{ij} 



23-d/2^{n+l)d/2 



k=l 



-Fal3\ 



5ii,i2,i3,i4 + Xl + X2 2 {n + l)d ' 



-F 





2 




- Xl - X2 




2 




+ Xl - X2 




2 


,22,23,14 


- Xl + X2 



2 {n + l)d 



where we introduced Zj^^ij.is,^ = ^ii + r^^ - r^g - r^^ and 2c|j^.} = ELi (i'i2;c-i " ^i^^ - Xfe 
Using 



dFai3{h,c^,a) 



-F^p{h,c^,a + 1), 



dR 



a/3 



dh, 



we may rewrite the answer as 



F{n) 

2n 22+n<i/2^(n+l)(i/2 



^ 1 n i^(xfe)dxfe 



X (r^, - - xij^ (ri3 - - X2j^ 
z. 



i^WiM + Xl + X2 2 + , 1 



V 


2 


^«l,*2,i3,i4 


Xl X2 




2 




+ Xl - X2 




2 


^il ,22 ,13 ,14 


- Xl + X2 



2 , 1 

Wi.a5 r 1 



Let us now derive the asymptotic form of at large distances for the case qo — Q. In this 
case -ff(x) decays fast at |x| > L and the expansion is derived in exactly the same way as 
the multi-pole expansion of the potential in electrostatics. We find 

, F(n) (/ H{^)x^d^f (/ i/(x)(ix)"-' 



2n 



22+n(i/2^2^(n+l)(i/2 



E 



X 



dbjdbs \ 



19 



note that 

J H{-K)x^di^ =-[{k + 2){k + d)D]-^ J x(x)a;^+2rfx . 
To complete the calculation we pass to study Fa/s- 



3. Computation of Faf3 



We now concentrate on F^/j defined by Eq. ()B12|) . By isotropy and incompressibility Fap 
satisfies 



= A' + 5' 



(B13) 



dhp " " ' ' b 
We consider the part of Fa/B that enters Fg^. Differentiating the above equation we find that 
for any two vectors v, w we have 



VaV.yWpWs — — — = — (V ■ W) 



db^dbs 



62 



(B' {d-2)B\ (B' 2B' 



(w ■ h) 



(v ■ b)^ 



fjB'V B' 85\(vb)^ 



/ ^.^2 (JB'\ {d-A)B' 2{d-A)B' 



X 



(v ■ w) (v • b) (w • b) 

62 



(B14) 



Thus our purpose is to find B. Taking the trace we find the additional equation that connects 
A and B: 

dA + B = DdT{a) [I{b) - I{b = 0)] , 

m^fdr, ^ 

J [(r - b)' + c2 

To perform angular integration conveniently we note that /(b) is a convolution of two 
spherically symmetric functions and its Fourier image is a product of spherically symmetric 
functions. To use Fourier transform techniques with logarithmic function we write Inr = 
-dr-yd^{^ = 0). Using 



dh 



exp [zb ■ q] _ (27r)'^/2 (cqf'"^ K,/^-a{cq) 



(62 + C2)" 



2"-iF(a)g 



d-2a 
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J a^-^r e 2) 



where K^, is a modified Bessel function of order 1/ — a — d/2, we find 



m 



d\. (27r)'='/2c-'^r((rf-0/2) 







\ '^"°6'^/2-i2-+«-ir(e/2)r(a) 

c2-+?r(d/2)r(a) V ^ 2 ' 2 ' 2 ' ' 

wliere Jd/2-1 is Bessel function of order d/2 — 1 and F is a liypergeometric function. Using 
I{b) and /(O) we find 

, , „ Ddr(iy)n'^/^ ( d\, 



d 

2' ^ 



where we used properties of hypergeometric function. One easily finds that B satisfies 

¥ db ~ 2^- \dbj[dij 

f(-.«44). (B15) 



The result of the differentiation over ^ can be written as a series using 



0, • =(^-^)'' 



so that 



We find for B 



DdT{u)n^/' - (_b_ 
2c2-(ci-l) ^^{d/2)k{d/2 + k) [ 6 



,2 



{d-l){d + 2)c^-^^ l^i/ + 1, 1, - + 2, --^ 



Next we find 



B' 4Dr(i/ + l)7r''/2 / d b 



2' 



(C^-I)(rf + 2)02*^+2 



FU+l,2,- + 2, 



-,2 
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It follows using 



that 



^ = — F{a + l,(3 + l,-f + l,z) 

oz 7 



h ) {d-l){d + 2){d + 4)02^^+2 \c\ 

d b^' 
2+^'-72 



xF(i/ + 2,3,- + 3,-- . (B17) 



Combining the results of the last two subsections one finds the answer for F^^ given in the 
main text. 
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